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PREFACE 


The  Air  Force  has  been  increasingly  concerned  with  logistics 
support  capability  under  very  demanding  and  dynamic  scenarios. 
Unfortunately,  many  measures  of  performance  used  in  evaluating 
logistics  support  and  methods  of  measurement  have  evolved  during 
relative  nond3mamic  peacetime  operations  and  conflicts  in  which 
demand  for  logistics  support  was  generally  stable.  These  "steady 
state"  methods  have  been  adapted  to  evaluating  dynamic  performance 
by  increasing  demand  for  support  to  higher  levels  and  assuming  that 
the  dynamic  scenarios  are  adequately  approximated  by  peacetime  steady 
state  demands  at  those  higher  levels.  This  assumption  could  lead  to 
the  overestimation  (and  in  some  cases,  underestimation)  of  wartime 
spares  requirements,  potentially  incorrect  war  reserve  spares  mixes 
for  squadrons,  and,  at  the  least,  inaccurate  capability  assessment  of 
alternative  logistics  policies. 

This  report  develops  techniques  of  transient  performance 
measurement  for  alternative  supply  and  maintenance  strategies.  A 
companion  document.  Part  II;  Component  and  Subcomponent  Spares 
Allocation  with  Dynamic  Models,  will  describe  the  use  of  these  models 
for  determining  spare  stockage  requirements. 

Rand-sponsored  research  led  to  the  dynamic  modeling  techniques 
described  in  this  report.  The  techniques  were  used  by  Carrillo, 
Hillestad  and  Lippiatt  in  "Maintenance  and  Supply  Performance 
Explorations,"  found  in  Appendix  E  in  [11],  and  were  further  expanded 
to  support  the  analysis  in  the  forthcoming  publication  "Assessing 
the  Capability  of  Planned  Reserve  Stocks  and  Spare  Engines  to  Support 
High  Sortie  Rates  in  a  Central  European  Contingency."  Because  the 
demands  for  the  Central  European  contlng'  ncy  are  expected  to  be  time- 
varying,  because  squadrons  are  collocated  with  potentially  shared 
resources,  and  because  the  availability  of  stockage  and  maintenance 
repair  resources  is  time-varying,  further  expansion  of  the  capability 
to  study  transitory  behavior  of  support  resources  was  necessary. 
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The  models  and  techniques  described  In  this  report  are  rather 
technical  and  a  background  In  probability  theory  Is  necessary  for  a 
full  understanding  of  the  necessary  assumptions,  development  of 
the  models,  and  some  Implications  of  their  use.  The  report  should 
be  of  greatest  Interest  to  those  specialists  Involved  In  studies  and 
analysis  of  dynamic  logistics  support  requirements,  those  Involved  In 
determining  wartime  stockage  (Including  spare  engines)  requirements, 
and  those  wanting  an  In-depth  understanding  of  the  models  used  by 
Rand  In  the  analysis. 

This  work  was  done  as  part  of  the  Project  AIR  FORCE  project 
"Strategies  to  Improve  Sortie  Production  In  a  Dynamic  Wartime  1 

Environment . "  I 
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I.  INTRODUCTION  AND  SUMMARY 


Much  of  the  past  work  regarding  recoverable  item  stockage  has 
concentrated  on  the  development  of  models  and  policies  for  systems  in 
steady  state.  The  well  known  work  of  Feeney  and  Sherbrooke  [3] 
regarding  an  (s-1,  s)  Inventory  policy  developed  distribution  functions 
for  items  in  service  and  measures  of  inventory  performance  based  on 
the  limiting  behavior  of  an  infinite  server  queue  expressed  by  Palm's 
theorem  [10].  Sherbrooke's  METRIC:  A  Multi-Echelon  Technique  for 
Recoverable  Item  Control  [14]  used  the  same  assumptions  regarding 
an  inventory  system  in  steady  state  with  constant  average  demand  rate 
and  service  rate.  More  recently  the  MOD-METRIC  model  of  Muckstadt 
[9]  expanded  the  capability  of  METRIC  by  considering  indentured 
components  but  also  assumed  that  recoverable  item  stockage  was  to  be 
provided  for  a  system  in  steady  state. 

These  models  and  variants  have  been  widely  applied  to  practical 
problems  of  inventory  management.  For  example,  the  Air  Force  uses 
METRIC-like  techniques  for  the  management  of  millions  of  dollars  of 
recoverable  spare  parts.  For  many  inventory  systems  including  the 
Air  Force  Supply  System  in  peacetime  the  assumption  of  steady  state 
behavior  is  both  convenient  and  adequate.  However,  there  are  important 
situations  in  which  the  transient  behavior  is  most  important.  A 
dramatic  example  of  this  is  the  potential  dynamic  behavior  exhibited 
by  demands  and  service  in  the  deployment  of  an  Air  Force  squadron  at 
the  onset  of  a  conflict.  Demands  for  components  may  suddenly  jump 
very  high  relative  to  the  previous  peacetime  activity  of  the  squadron 
and  then  decay  gradually  or  abruptly  due  to  the  attrition  of  the 
aircraft  in  the  squadron.  Meanwhile  the  initial  service  rate  may 
be  near  zero,  as  the  already  deployed  unit  awaits  the  airlift  of  the 
specialized  personnel  and  test  equipment  to  repair  broken  components, 
gradually  Increase  to  its  full  wartime  service  capability,  and 
suddenly  drop  to  near  zero  again  as  a  result  of  damage  in  an  airbase 
attack. 
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Two  problems  regarding  recoverable  Item  inventory  faced  by  a 
decisionmaker  for  the  nonstationary  case  are  the  same  as  those  faced 
in  the  steady  state  case,  that  is,  to  determine  how  much  spare 
stock  to  provide  against  stockouts  and  to  determine  what  level  of 
performance  can  be  achieved  with  a  given  level  of  investment  in  spare 
stock.  Clearly,  the  limiting  behavior  of  the  dynamic  case  is  not  of 
interest  when  that  limit  is  reached  at  times  far  beyond  the  period  of 
interest.  Performance  measured  with  models  which  only  consider 
limiting  behavior  will  tend  to  be  inaccurate  and  stockage  provided 
against  limiting  behavior  may  be  either  inadequate  or  excessive. 

Gaver,  Lehorsky,  and  Perlas  [5]  have  described  certain  service 
systems  with  transitory  demands.  Miller  [11]  has  described  a  real¬ 
time  management  approach  for  the  distribution  of  stockage  among 
locations  affected  by  time  varying  demands.  The  primary  purpose 
of  this  Note  (Part  I)  is  to  provide  an  Integrated  approach  to  inventory 
performance  measurement  for  a  given  stockage  allocation  for  systems 
with  nonstationary  demands  and  service  rates.  Part  II  will  describe 
certain  stock  allocation  algorithms  for  the  nonstationary  case. 

A  secondary  purpose  of  the  Note  is  to  suggest  approaches  to 
certain  aspects  of  recoverable  item  repair  and  supply  which  currently 
cause  significant  deviation  between  practice  and  theory,  even  in  the 
case  of  stationary  demands  and  service  rates.  In  particular,  we  will 
address  the  problems  of  cannibalization  of  primary  recoverable  items 
by  providing  models  for  different  degrees  of  cannibalization. 

Cannibalization  of  components  is  common  practice  in  many  service 
systems,  particularly  in  cases  where  the  equipment  to  be  repaired  is 
highly  modular  and  interchangeable  with  a  minimum  of  effort.  In 
effect  it  provides  an  additional  source  of  spare  components  when 
stock  is  low  or  service  times  are  long.  For  example,  consider  the 
case  of  an  aircraft  not  operational  due  to  the  shortage  of  one  item. 
While  that  aircraft  is  nonoperational  it  also  provides  a  hedge  against 
other  aircraft  becoming  nonoperational  due  to  shortages  of  other  items. 
In  the  same  fashion  Indentured  items  awaiting  service  provide  a 
secondary  source  of  supply  of  subcomponents.  The  METRIC  and  MOD-METRIC 
techniques  do  not  consider  the  effect  that  cannibalization  has  on 
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inventory  system  performance  and  therefore  tend  to  overestimate 
stockage  for  systems  which  permit  cannibalization.  We  will  describe 
how  cannibalization  can  be  considered  in  evaluating  performance. 

We  will  also  describe  a  group  of  performance  measures  which 

attempt  to  predict  the  effect  that  inventory  and  service  policy  has  on 

organizational  performance.  Some  of  these  measures  have  been  developed 

for  the  special  case  of  aircraft  support  in  which  the  organizational 

goal  is  to  provide  enough  operational  aircraft  to  perform  missions. 

(However,  it  should  be  quite  easy  to  adapt  them  to  support  of  other 

equipment.)  Because  the  measures  deal  with  end  item  performance 

* 

and  demands  they  are  scenario  dependent.  For  example,  stockout 
of  aircraft  components  may  not  be  important  if  the  aircraft  are  not 
required  to  fly  large  numbers  of  flights  since  only  a  few  operational 
aircraft  are  needed  to  perform  the  activity.  Thus,  we  will  provide 
a  model  with  the  potential  to  evaluate  inventory  policy  in  the 
context  of  organizational  goals. 

Finally,  in  Appendix  B  we  have  indicated  certain  shortcomings 
in  our  models  and  other  currently  used  models  which  have  not  been 
overcome  as  far  as  we  know.  We  have  also  suggested  the  possible 
effects  of  these  shortcomings. 

Certain  parts  of  the  models  in  this  paper  have  been  used  for 
evaluation  of  ship  and  shore  support  of  the  carrier  aircraft  in 
the  Defense  Resource  Management  Study  by  Rice  [11],  for  analysis 
of  centralized  maintenance  support  of  the  A-10  aircraft  in  the 
same  study,  and  for  analysis  of  recoverable  item  stockage  requirements 
for  various  levels  of  potential  wartime  flying  demands  of  the  F-4 
aircraft. 

The  next  section  gives  the  derivation  of  a  basic  results 
regarding  the  time  dependency  of  recoverable  items  in  a  queue.  This 


See  Appendix  A  for  the  methodology  needed  to  include  scenario 
information  in  the  model. 
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is  followed  by  a  section  describing  a  group  of  time  dependent  measures 
of  performance. 
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II.  THE  BASIC  RESULT  -  TIME  DEPENDENCY  OF  NUMBER  OF 
COMPONENTS  IN  A  QUEUE 


In  this  section  we  will  obtain  the  probability  distribution  of 
the  number  of  recoverable  items  in  service  when  the  demand  rate  for 
service  and  the  service  rate  are  nonstationary  and  there  is  sufficient 
reserve  service  capacity-  This  distribution  will  then  be  used  in 
calculating  time  varying  performance  measures  during  periods  of 
transitory  demand  and  service.  First,  we  review  a  related  steady 
state  result. 

The  well  known  queuing  theorem  due  to  Palm  ([10],  see  also  [15]) 
states  that  when  demands  arrive  according  to  a  Poisson  process  with  rate 
X  and  the  service  times  have  an  arbitrary  probability  distribution 
Independent  of  the  demand  process  with  mean  T  then  the  steady  state 
probability  of  k  units  in  service  under  the  assumption  of  slack  repair 
capacity  is  given  by  a  Poisson  distribution  with  parameter  XT,  that  is. 


probability  of  k  units  in  service  =  Pj^ 


(XT)V^" 

k! 


Palm  provides  a  similar  result  for  the  lost  sales  case  in  which  demands 
stop  when  there  is  stockout.  Feeney  and  Sherbrooke  [4]  extend  these 
results  for  the  case  of  the  compound  Poisson  distribution.  All  of 
these  results  are  valid  only  in  the  limit  and  when  both  the  demand 
and  service  distribution  are  stationary.  We  will  derive  a  similar 
result  for  time  varying  service  and  demand  process. 

We  begin  our  derivation  for  the  nonstationary  case  with  the 
nonhomogeneous  Poisson  process.  A  demand  process  is  said  to  be  a 
nonhomogeneous  Poisson  process  with  intensity  function  m(t) ,  t  ^  0 
if 


i.  The  number  of  demands  existing  at  time  t  =  0  is  zero, 
ii.  The  numbers  of  demands  in  disjoint  time  increments 
are  independent  of  each  other. 
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iii.  The  probability  of  more  than  one  demand  in  an  increment 

becomes  infinitesimally  small  as  the  increment  gets  small. 
(A  compound  Poisson  process  relaxes  this  assumption.) 
iv.  The  probability  of  one  demand  in  any  increment  is  given 
by  the  intensity  function  m(t)  times  the  length  of  the 
increment  as  the  increment  gets  small. 

Then,  if  we  let 


N(t)  =  Number  of  demands  by  time,  t 


Y(t) 


m(s)ds 


it  can  be  shown  (see  Ross  [13],  for  example)  that 

Probability  of  k  demands  in  the  interval  t,  t  +  s  is 

P[S(t  +  =)  -  N(t)  -  kl  - 
for  k  ^  0 

Y(t)  is  called  the  mean  value  function  of  the  process. 

The  following  result  is  stated  in  an  exercise  in  [12]. 

Theorem:  Given  N(t)  =  k,  the  unordered  set  of  arrival  times  in  a 
nonhomogeneous  Poisson  process  with  mean  value  function  Y(t)  have  the 
same  distribution  as  k  independent  and  identically  distributed  random 
variables  having  distribution  function 

X  £  t 
X  >  t 


i 
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Theorem.  (Nonhomogeneous  Poisson  Queue) 

Let  X(t)  denote  the  number  of  customers  in  the  system  at  time  t 
when  arrivals  occur  according  to  a  nonhomogeneous  Poisson  process 
with  intensity  function  m(t)  and  mean  value  function 


■Y(t) 


o 


Let  the  nonstationary  service  distribution  be  F(s,t)  (the  probability 
that  a  service  started  at  time  s  will  be  completed  by  time  t) .  Then, 
with  the  additional  assumptions  that  the  repair  process  is  independent 
of  the  arrival  process  and  slack  repair  capacity,  X(t)  has  a  Poisson 
distribution  with  mean 


F(s,t)dY(s) 


F(s,t)m(8)ds 


where 


F(s,t)  =  1  -  F(s,t). 


Proof : 

Conditioning  on  N(t)  we  obtain  the  probability  that  X(t)  =  j. 


P(X(t)  =  j) 


P(X(t)  =  jlN(t)  =  k)  -  ^ 

k=0 


A  version  of  this  theorem  as  well  as  a  sketch  of  a  proof 
(different  from  the  one  given  here)  was  first  seen  by  the  authors 
in  some  unpublished  notes  by  Gordon  Crawford.  Also,  see  Crawford  [2] 
and  Gaver  [6]. 
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If  there  have  been  k  arrivals  by  time  t,  then  the  probability 
that  an  arbitrary  one  of  these  customers  will  still  be  present  at  t 
is  given  by 


t_ 

i-  (s,t)d 


This  is  obtained  by  conditioning  on  the  time  of  arrival  of  an  arbitrary 
customer  and  then  unconditioning  using  the  distribution  G(x)  given 
earlier. 

Now, 


(j)  j  =  0.  1,  ....  k 


P(X(t)  =  j|N(t))=  k  = 


j  >  k 


Uncondltionlng  on  N(t),  we  get 


p(x(t)  =  j)  =  P^d-P)’"'^ 

k=j  ' 


-3 


k! 


=  e"^ 


OO 

(t)  (p*Y(t))^  \  ’’  [  (l-p)Y(t)]^~^ 

j!  (k-j)! 


k=j 


^  g-Y(t)  (p»Y(t))^  ^(l-p)Y(t) 
_  „-PY(t)  [p»Y(t)]^ 

J! 


Furthermore,  if  we  let  X(t)  and  Y(t)  be  the  number  of  customer <  Jn 
service  in  two  independent  processes  it  can  be  shown  that  the  total 
number  in  service  is  given  by  the  Poisson  distribution. 


P(X(t)  +  Y(t)  =  k)  = 


[X^(t)+XY(t)]’' 


-[Xj^(t)+XY(t)] 

e 


k! 


We  will  use  the  following  example  to  give  some  perspective  on  these 
results. 


AN  AIRCRAFT  ENGINE  EXAMPLE 

Assume  that  aircraft  engines 
average  rate  m(t)  where 


at  an  airbase  demand  service  at  an 

t  <  T 

t  >  T 


and  that  t  reflects  a  change  in  the  base  flying  program  causing  a 
change  in  the  rate  of  demand  of  service.  Assume  that  an  arbitrary 
service  demand  has  a  probability  f  of  being  satisfied  at  the  base  and 
must  be  satisfied  at  a  remote  location  with  probability  1-f.  Assume 
further  that  the  service  time  distribution  functions  for  base  and 
remote  locations  are  exponential  with  average  rates  p  and  p  .  Finally, 
assume  that  engines  demanding  service  at  the  remote  location  cannot 
leave  the  base  for  the  first  T  (assume  T  <  x)  days  of  the  period  we 
are  considering.  Then  we  have,  for  base  service 
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and  for  remote  service 


s  <  t  <  T 


Fj^(s,t)  = 


-Pj(t-s) 


8  <  T  <_  t 
T  <  s  <  t 


The  average  number  of  engines  In  the  base  repair  process  Is  then 
given  by 


r 


Xg(t)  =  I  Fg(s,t)m(s)ds 


/t  -M  (t-s)  m  -y  t 

e  m^fds  =  f  •  —  (1-e  ) 

B 
o 


/t  -y  (t-s) 

e  m^^fds  +  / 


-y  (t-s) 
e  m^fds 


for  t  <  T 


for  t  >  T 


f  •  —  (1-e  )  e 

y 


ra  -y  (t-T) 

+  f  •  ^  (1-e  ®  ) 


for  t  >  T 


and  the  average  number  of  engines  In  repair  at  the  remote  location  or 
awaiting  repair  at  the  remote  location  Is 


XR(t) 


■r 


li 


Fj^(s,t)m(s)ds 


'-rV 


=  J"  1  •  inj^(l-f)ds  =  (l-f)mj^t 


for  t  <  T 


ri  -y  (t-T)  /t  -y  (t-s) 

X^(t)  =1  e  m^(l-f)ds  +  /  e  m^(l-f)ds 


-y  (t-T)  m  -y  (t-T) 

=  (l-f)m  T  e  +  (1-f)  —  (1-e  )  for  T  <  t  <  t 

/•T  -y  (t-T)  ri  -y  (t-s) 

XR(t)  =  /  e  inj^(l-f)ds  +  /  e  m^(l-f)ds 


/•t  -y  (t-s) 

+  I  e  m2(l-f)ds 


=  (l-f)mj^T  e 


-yj^(t-T) 


m  -y  (t-T)  -y  (t-i) 

+  (1-f)  —  (1-e  ^  )  e 

^R 


+  (1-f)  —  (l-e~^R^*^~'^^) 
^R 


for  T  <  t 


The  average  total  engines  In  service  is  then  given  by 


A„(t)  =  X  (t)  +  A_(t) 
1  B 


and  the  distribution  of  engines  in  or  awaiting  service  as  a  function 
of  time  is 


P(Xg(t)  +  Xj^(t)  »  k)  = 


k 

A^(t)''  e  ^ 


t 

r 
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Figure  1  shows  the  mean  number  being  serviced  as  a  function  of  time. 
Note  that  as  t  gets  larger  this  result  converges  to  the  steady  state 
result  predicted  by  Palm's  theorem, 


m.  m* 

't<")  ■  *  4  ^  4' 


Note  also  that  the  dynamic  result  gives  considerably  more  information 
regarding  the  distribution  of  engines  in  service  during  the  period  of 
transient  behavior.  It  can  be  used,  for  example,  to  determine  the 
sensitivity  to  base  and  remote  repair  rates,  the  effect  of  transition 
to  a  different  flying  rate  (with  subsequently  different  demand  rate) 
and  the  effect  of  various  delays  in  shipping  engines  off  base. 
Furthermore,  it  can  be  used  to  determine  the  time  dependency  of  the 
probability  of  stockout  based  on  a  time  varying  supply  of  spare 
engines. 

In  this  example  we  assumed  that  there  were  no  engines  requiring 

ic 

service  Initially.  However,  by  superposition  it  is  possible  to  show 
that  if  there  was  an  average  of  engines  requiring  service 
initially,  the  mean  number  in  service  at  time  t  is 


A(t) 


F(s,t)m(s)ds  +  F(0,t)rQ 


provided  that  service  times  have  an  exponential  distribution. 

Appendix  A  gives  a  number  of  closed  form  solutions  for  X(t) 
given  various  assumptions  about  F(s,t)  and  m(s) . 


Consider  two  separate  repair  processes,  one  which  operates  prior 
to  t  *  0,  and  another  which  operates  after  t  »  0.  The  second  term 
in  the  equation  above  represents  the  servicing  of  all  demands  which 
have  occurred  and  not  been  serviced  by  time  t  =  0. 


-Average  number  of  engines  being  serviced  as 
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Once  an  expression  or  value  of  X(t)  is  determined  It  Is  possible 
to  compute  other  measures  of  performance  such  as  the  probability  of 
stockout  given  a  spare  stock  level,  S(t),  at  time  t.  The  next  section 
derives  a  number  of  performance  measures  for  the  nonstationary 
inventory  system. 
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III.  TIME  DEPENDENT  MEASURES  OF  PERFORMANCE 


In  this  section  we  describe  a  number  of  performance  measures  of 

* 

a  nonstationary  inventory  system  ,  some  of  which  are  the  well  known 
measures  such  as  average  backorders,  while  others  are  more  specifi¬ 
cally  related  to  system  supporting  activity  of  vehicles  such  as 
aircraft.  Let  X^(t)  represent  the  mean  number  of  a  single  recoverable 
item,  1,  undergoing  service  as  given  in  the  previous  section  and  let 
S^(t)  be  the  amount  of  spare  stock  provided  for  that  item  at  time  t. 
(We  assume  here  that  S^(t)  is  a  deterministic  quantity.)  Then 
several  Inventory  system  measures  are: 

a)  Probability  of  stockout  at  time  t,  PO^(t) 

This  is  the  probability  that  demands  exceed  supply 
and  is  given  by 


PO^(t)  = 


-Xj(t) 


k=S^(t)+l 

S^(t) 


=  1 


-x,(t) 


(t) 


k  e 


k=0 


b)  Fill  rate  at  time  t,  FR^(t) 

This  is  the  probability  that  a  demand  that  has  occurred 
can  be  filled 


FR^(t) 


S^(t)-1 

E 


k=0 


X^(t) 


-\(t) 

k  e _ 

k! 


See  [1]  for  steady  state  supply  performance  measures. 


c)  Average  backorders  at  time  t,  EB^(t) 

This  gives  the  average  level  of  shortages  of  item  i 


at  time  t 


J^(t)  =  ^  [k  -  S^(t)}  -i 


A,(t)’"e  ^ 


k=S^(t)+l 


S.(t) 


Ea  , 

(S^(t)  -  k)  — 


k 

A,(t)‘"e  ^ 


d)  Variance  in  backorders  at  time  t 


\  T  A.(t)  e  „ 

VBi(t)  =  >  [k  -  S^(t)]^  - - EB^(t)^ 

k=^t)+l 


=  A.(t)  +  [A^(t)  -  S.(t)]- 


S^(t) 


k 

A,(t)‘'e 


-  EB^(t)‘ 


e)  Average  nvimber  of  systems  NMCS  (Not  Mission  Capable  for 
Supply  reasons)  without  cannibalization,  EN(t) 

This  measure  gives  the  average  combined  effect  of  recoverable 
item  shortages  on  the  major  system  (such  as  aircraft)  those 
items  support.  Here  we  assume  that  NA(t)  represents  the 
number  of  such  major  systems  which  are  supported  at  time  t. 

We  further  assume  that  there  are  N  types  of  recoverable 
items  required  on  each  system  and  that  the  shortages  of 
any  one  of  these  items  will  make  the  system  nonoperatlonal. 

In  the  noncannlballzatlon  case  we  assume  that  shortages  of 
components  cannot  be  consolidated  among  the  systems. 
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The  probability  that  an  arbitrary  system  has  a  shortage  of 
item  i  when  there  are  k  shortages  of  i  is 

k 

NA(t) 

and  therefore,  the  probability  that  an  arbitrary  system  has 
a  shortage  of  item  i  is 


EB^(t) 

NA(t) 

The  probability  that  a  system  is  down  due  to  shortages  of 
some  item  is 


1 


EB^(t) 

NA("t) 


and,  finally,  the  expected  number  of  systems  not  operational 
at  time  t  is 


EN(t)  =  NA(t) 


EB^(t) ^ 
NA(t)  / 


This  derivation  required  the  assumption  that  only  one 
component  of  given  type  was  on  each  system.  If  this  is  not 
the  case  let  be  the  quantity  of  item  i  per  system.  If 
no  consolidation  of  shortages  occurs  then 


EN(t)  =  NA(t)  I  1 


-n 


Qj^NA(t)  -  y 


PB^(y) 


Q^NA(t) 


where  PB^(y)  is  the  probability  that  item  i  has  y  shortages 
at  time  t 


PB^(y)  = 


Si(t)  w 

e  X^(t) 


-X  (t)  S  (t)+y 
e  (t) 


(S^(t)+y)! 


y  =  0 


y  >  0 


f)  NORS  with  complete  cannibalization 

Here  we  assume  that  shortages  of  all  components  are 
consolidated  to  make  the  smallest  number  of  nonoperational 
systems.  Let  P^(j)  be  the  probability  that  shortages 
of  the  i*"^  item  are  less  than  or  equal  to  j .  Then 


S.(t)+j 


P^j)  =  53 


k 

X^Ct)*"  e  ^ 


Let  P(j)  be  the  probability  that  the  number  of  nonoperational 
(NORS)  systems  is  less  than  or  equal  to  j .  Then 


N 

P(j)  = 
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If  there  is  more  than  one  item  of  a  type  on  each  system 
we  again  employ  as  the  number  of  item  i  per  system  and 
obtain 

N 

pa)  =  I””!  P^CQ^-j)- 

i=l 


The  average  number  of  NORS  systems  with  full  cannibalization 
is  then 


NA(t)-l 

ENc(t)  =  2  "  pa)]  • 

j-0 

The  NORS  distribution  function  with  full  cannibalization 
is 


PNj(t)  =  P(j)  -  pa  -  1) 

which  then  can  be  used  to  give  the  variance  in  the  number 
of  NORS  systems. 


VN(t)  = 


NA(t) 

Z 


j=l 


j^PNj(t) 


-  EN^Ct) 


g)  NORS  with  partial  cannibalization 

Here  we  assume  that  some  items  are  relatively  easy  to 
cannibalize  and  that  some  are  so  difficult  to  remove  or  Install 
that  it  is  not  desirable  to  cannibalize  them.  Let 


I 

c 


I 

n 


item  1  is  cannibal izable| 
item  i  is  not  canniballzable 
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We  first  compute  the  probability  of  exactly  k  NORS  systems 

due  to  items  in  the  set  1  .  This  is  just  the  PN.(t)  derived 

1 

earlier  using  P(j)  and  where  P(j)  is  computed  only  using  the 
i  in  I^.  We  will  denote  this  probability  PN?(t).  The 
probability  that  an  arbitrary  system  is  operational  considering 
only  shortages  of  cannibalizable  items  is  then 


NA(t) 


E 

4=1 


-  j  ■- 

NA(t) 


PNj(t) 


EN^(t) 

NA(t)‘ 


The  probability  that  an  arbitrary  system  is  operational  after 
shortages  of  noncannibalizable  components  only  is 


n 


assuming  Q^,  the  quantity  of  item  i  per  system,  is  one  for  each 

part,  i,  in  the  set  1^.  The  probability  that  an  arbitrary  system 

is  not  operational,  assuming  independence  of  demands  and  that  the 

cannibalization  of  the  items  in  I  takes  place  with  no 

c 

information  of  the  failed  items  belonging  to  I^,  is 


1  - 


1 


EN^(t) 

NA(t) 


EB^(t) 

NA(t)^ 


and  the  expected  number  of  nonopera tional  systems  given 
partial  cannibalization  and  given  NA(t)  total  systems  is 


ENP(t) 


NA(t) 


EN^(t) 

NA(t) 


n 

iel 


EBjt) 


J 
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EN^(t)  Is  an  overestimation  of  the  expected  NORS  if  the 
information  describing  which  noncanniballzable  items  had  failed 
is  used  (since,  in  this  case,  the  strategy  would  be  to  move 
"holes"  due  to  cannibal Izable  items  to  those  aircraft  with  holes 
due  to  noncanniballzable  items.  We  did  not  do  this  in  our 
derivation) . 

In  the  NORS  expressions  found  in  sections  e,  f,  and  g,  the 
reader  may  observe  that  backorders  and  shortages  probabilities 
have  been  computed  as  though  demands  were  Independent  of  the 
number  of  systems  nonoperational.  In  most  systems  the  demand 
rate  lowers  as  the  number  of  systems  nonoperational  increases 
due  to  lower  activity  of  the  systems.  Unfortunately,  the 
independent  increment  assumption  of  the  Poisson  process  is  then 
violated.  Usually,  with  spare  stock  provided  for  the  appropriate 
level  of  activity,  the  percentage  of  systems  nonoperational  is 
low  so  that  there  is  a  very  low  probability  that  failures  exceed 
spare  stock  plus  the  number  of  systems.  For  cases  in  which  NORS 
percentages  get  high,  it  is  possible  to  use  the  lost  sales 
distribution  with  the  demand  turned  away  when  the  number  of  demands 
exceeds  spare  stock  plus  the  number  of  systems, 
h)  Probability  of  meeting  system  demands  and  number  of  system 

demands  met 

This  group  of  performance  measures  is  for  the  situation  in 
which  it  is  desirable  to  have  enough  systems  operational  to 
perform  a  certain  demanded  level  of  activity.  In  the  case  of 
aircraft  systems  this  might  be  the  number  of  operational  airframes 
which  permits  a  desired  level  of  flying  activity.  Let  B(t)  be 
the  maximum  number  of  allowable  NORS  systems  which  still  permits 
the  system  demands  to  be  met  and  let  D(t)  be  the  level  of  system 
demands  at  time  t.  Let  r(t)  be  the  amount  of  system  demand  which 
can  be  satisfied  by  a  single  system  at  time  t.  In  the  aircraft 
case  D(t)  would  be  the  total  desired  sorties  per  unit  time  while 
r(t)  would  be  the  maximum  number  of  sorties  per  unit  time 
achievable  by  a  single  aircraft. 


1 
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Note  that 


B(t)  =  NA(t)  - 


D(t) 

r(t) 


where  [x]  is  the  ceiling  of  x,  that  is  the  smallest  integer 
y  so  that  y  ^  x,  and  D(t)  ^  r(t).  NA(t)  is  assumed. 

■t.  The  probability  of  meeting  system  demands  is  then 


PD(t)  =  P(B(t)) 


That  is,  the  probability  that  system  demands  are 
met  is  just  the  probability  that  the  number  of 
nonoperational  systems  is  less  than  or  equal  to 
B(t). 

ii.  The  expected  number  of  demands  met  given  k 
nonoperational  system  is 

I  D(t)  if  k  <  B(t) 

ES(t)/k  =  < 

/  r(t)(NA(t)-  k)  if  k  >  B(t) 
Unconditioning 

NA(t) 

ES(t)  =  D(t)P(B(t))  +  y"'  r(t)(NA(t)  -  k)PN  (t) 

k=B(0+l 

where  P(j)  and  PNj^(t)  were  described  earlier. 
iii.  The  distribution  of  the  number  of  demands  met 
is  given  by 

if  k  *  D(t) 

if  k  -  r(t)(NA(t)  -  j); 
j  -  B(t)  +1,  ....  NA(t) 

otherwise 


PS^(t) 


P(B(t)) 


PNj(t) 
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iv.  The  variance  In  demands  met  is  given  by 

NA(t) 

VS(t)  =  D^(t)  P(B(t))  +  r^(t)  (NA(t)  -  k)^  PN  (t)  -  ES^(t) 

k-W+l 

Note  that  the  latter  group  of  measures  bring  into  play, 
in  a  deteministic  sense,  additional  time  varying  parameters 
such  as  the  number  of  systems,  NA(t),  the  demands  on  all  systems, 

D(t),  and  the  portion  of  demand  a  single  system  can  satisfy, 
r(t).  In  the  deployment  of  an  aircraft  squadron,  NA(t)  might 
represent  a  time  phased  deployment  of  the  aircraft  as  well 
as  attrition  due  to  losses  in  a  conflict.  The  demand  D(t) 
would  be  the  time  dependent  demand  for  sorties  as  an  engagement 
proceeds  and  r(t)  would  be  time  varying  because  of  a  short¬ 
term  capability  to  provide  surges  of  flying  activity. 
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APPENDIX  A 

NONSTATIONARY  SERVICE  DISTRIBUTIONS 


In  this  section,  we  will  present  various  non-stationary  service 
distributions  which  are  of  interest  in  evaluating  the  sensitivity 
of  a  recoverable  item  inventory  system  where  the  demands  for  items 
arrive  according  to  a  nonhomogeneous  Poisson  process,  implying 
changing  rates  of  demand  for  service.  Under  the  assumption  of 
sufficient  slack  service  capacity  and  the  independence  of  the  demand 
and  service  processes,  it  was  shown,  in  the  main  sections  of  this 
report,  that  for  such  systems  the  number  of  items  in  service  at  time 
t  is  given  by  the  integral 


A(t)  =  J  F(s, t)m(s)ds 

o 

where  F(s,t)  is  the  probability  that  a  service  started  at  time  s 
will  not  be  completed  before  time  t  and  m(s)  is  the  Intensity  function 
of  the  Poisson  process. 

The  service  distributions  to  be  studied  have  one  main  property 
in  common:  the  service  is  either  nonexistent  for  all  or  some  period 
of  time  or  a  service  process  is  available  with  service  times  that 
have  an  exponential  distribution  with  rate  w.  This  property  will 
allow  us  to  write  X(t)  in  terms  of  two  functions: 


Y(t) 


m(s)ds. 


the  mean  value  function  at  time  t  of  the  nonhomogeneous  Poisson 
process,  and 


,(t) 


/ 


-Vl(t-8) 


m(s)ds 


J 
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whlch  is  discussed  next. 


A)  Basic  Case: 


Al.  Uninterrupted  service  process  with  exponential  service 

times,  l.e.,  F(s,t)  =  1  -  e  for  0  £  s  £  t. 

A2.  Constant  demand  rate  m,  that  Is 


m(s)  = 


s  <  0 


s  >  0 


Clearly,  In  this  case 


-r 


m(s)ds  =  m  t. 


The  number  of  items  in  service  is  a  Poisson  random 
variable  with  mean 


A(t)  =  /  e'^^’''®^m{s)ds  =  -K  (t) 

Jo  ^ 


For  the  assumed  m(t),  wc  get 


.^(t)  =  H 


Note  that  for  large  t,  when  the  system  has  been  operating 
for  a  long  time,  we  have,  that  the  number  of  Items  In 
service  approximates  — ,  the  steady  state  result. 


B)  A  rfore  General  Demand  Rate  for  Service: 
Bl.  Same  as  Al. 
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Piecewise  constant  m(t),  l.e.. 


m(s)  = 


s  <  0 


i-1  <  t  _<  i,  for  i*l,  t,  t  assumed  Integral 


so  that  the  mean  value  function 

y(t)  =  f  m(s)ds  =  m^  for  integral  t  >  0. 

1=1 

Here,  the  mean  number  in  service,  A(t),  becomes 


A(t)  = 


for  t  =  2,  we  have 


-y(2-s)  , 

e  m2ds 


in.  lu. 

=  (l-e'^)e"^  +  77^  (1-e  ") 


m„ 


For  Integral  t  ^  1 


X(t)  *  $y(t) 


~2^  V 


-u(t-i) 


i-1 


J 
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C)  Service  Process  Unavailable  for  Some  Initial  Time  Period: 

Cl.  No  service  available  on  [0,  x],  i.e. ,  up  to  time  x; 

thereafter,  have  exponential  service  times  with  mean 
C2.  General  intensity  function  m(s)  for  0  s  ^  t. 

In  this  case, 

^  Y(t)  for  t  £  X 

Ht)  =  I 

|Y(x)e"^^‘^“'^^  +  $^(t)  -  *J>^(t)  for  t  >  X 


D)  Subsequent  Loss  of  Repair  Process: 

Dl.  Exponential  service  time  (rate  p)  on  [0,  x] ,  i.e.,  before 
time  T,  after  which  the  system  completely  loses  its 
service  capability. 

D2.  Same  as  C2. 

Here 

(t)  for  t  <  X 

$^(x)  f  Y(t)  -  y(x)  for  t  >  X 

E)  Service  Times  with  a  Fixed  Time  Component: 

El.  Service  times  made  up  of  a  fixed  component  d  plus  an 
exponential  random  component  with  mean  — . 

E2.  Same  as  C2. 

The  result  is 


A(t) 


Y(t)  t  £  d 

'{'^(t-d)  +  Y(t)  -  Y(t-d)  d  <  t 


F)  Service  Times  with  Fixed  Ccmiponent  Plus  Loss  of  Service : 

FI.  Same  as  El  except  chat  after  time  x  the  system  permanently 
loses  its  service  capability. 
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F2.  Same  as  E2. 

Get 

t  £  d 
d  <  t  £  T 
T  <  t;  t  >  d 


!Y(t) 

$^(t-d)  +  Y(t)  -  Y(t-d) 

<3>^(max(0,  T-d))  +  Y(t)  -  Y(max(0,  T-d)) 


G)  Discrete  Change  in  Service  Rate: 

Gl.  Assume  exponential  service  times  with  rate  on  [0,  » 

i.e.,  prior  to  time  t^^,  after  which  the  rate  permanently 
changes  to  y^. 

G2.  Same  as  C2. 

We  obtain 


(t) 


A(t) 


-y„(t-tp  __  -y,(t-t  )  _  _ 

'  (t)  -  e  (t,)  t,  <  t 


H)  Two  Discrete  Changes  in  Service  Rate 

HI.  Assume  exponential  service  times  with  rates  y2»  y2» 

in  the  respective  time  intervals  [0,  t^] ,  (t^,  t2]»  ^2’ 
where  0  <  t^  <  t^- 
H2.  Same  as  C2. 

The  result  is 


A(t) 


(t)  t  <  t 

^  ^ 
-y  (t-t  )  _  -y  (t-t  ) 

e  $  (t.)  +  (t)-e  •  $  (t-) 

1  2  ^2 
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e  e  $  (t  ) 

\  ^ 


+  e 


-y3(t-t2) 


$  (t.)  -  e  $  (t.) 

^2  2  ’^2  ^ 


*^2  ^  ^ 


-vi_(t-t  )  _ 

+  $  (t)  -  e  ^  ^  $  (t  ) 

M3  li3  ^ 


I)  Two  Discrete  Changes  in  Service  Rate  and  Service  Times  with 
Fixed  Component 

11.  Same  as  HI  but  in  addition,  there  is  a  delay  of  d  units 
of  times  prior  to  an  item  starting  its  service. 

12.  Same  as  C2. 


A(t)  = 


Y(t) 


t  <  d 


$  (t-d)  +  Y(t)  -  Y(t-d) 


d  <  t  < 


-M  (t-t  )  _ 

e  (max(0,  tj^-d)) 


+  $  (t-d)  -  e 


-y^Ct-tj^) 


(max(0,  t^-d))^for  <  t  <  t2;  t 


+  Y(t)  -  Y(t-d) 


-y3(t-t2)  _ 

e  e  (max(0,  t^^-d)) 


+  e 


-P3(t-t2)| 


I*  (max(0,  t„-d))  -  e  $  (max(0,  t,-d))| 


+  $  (t-d)  -  e 


-P3(t-t2) 


'^jj3(niax(0,  t2-d))  +  Y(t)  -  Y(t-d) 


for  t-  <  t;  t  >  d 


As  it  becomes  clear,  further  assumptions  would  Imply  more 
complicated  expressions  which  may  lead  to  numeric  integration.  However, 
the  various  service  distributions  presented  allow  a  quite  general 
sensitivity  testing  of  an  inventory  system. 
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APPENDIX  B 


SOME  UNRESOLVED  PROBLEMS 


The  development  of  models  and  techniques  for  transitory  behavior 
in  the  body  of  this  Note  required  certain  assumptions  that 
are  sometimes  violated  to  the  extent  that  the  models  are  no  longer 
good  approximations  of  the  situation  modeled.  The  most  important 
assumptions  are: 

1.  Sufficient  slack  service  capacity 

2.  Independence  of  the  service  and  demand  process 

3.  Poisson  arrivals 

The  first  assumption  implies  that  each  recoverable  item  demanding 
service  immediately  receives  service  with  average  service  time  based 
on  the  fuiiction  F(s,t).  In  general,  the  service  time  is  dependent 
on  the  number  of  servers,  the  scheduling  of  repair,  and  the  number  of 
components  already  placing  demands  on  servers.  With  respect  to 
aircraft  recoverable  components  the  number  of  servers  may  be  limited 
by  the  availability  of  specialized  test  equipment  or  specialized 
personnel.  Furthermore,  the  test  equipment  and  personnel  are  generally 
used  to  service  several  or  many  different  kinds  of  recoverable  items. 
Queuing  results  for  finite  numbers  of  servers  are  available  but  have 
the  drawback  that  they  do  not  account  for  the  cross  training  and 
sharing  of  servers  for  different  types  of  repair.  Furthermore,  they 
do  not  lead  to  the  Poisson  result  described  in  this  Note 
and  consequently  do  not  allow  the  separation  and  combination  of 
different  repair  processes  or  "pipeline." 

The  "steady  state"  stockage  calculations  in  use  by  the  Air  Force 
today  use  an  average  resupply  time  as  the  average  service  time.  This 
average  resupply  time  Includes  waiting  time  in  service  queues  and  is 
a  historical  average  of  the  length  of  time  elapsed  between  when  it 
is  checked  in  for  repair  and  checked  out  of  repair.  An  infinite  server 
model  using  this  time  as  the  average  service  time  is  then  used  for 
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recoverable  Item  stockage.  This  approximation  is  reasonable  as  long 
as  the  system  remains  strictly  in  steady  state  or  when  there  is 
enough  slack  capacity  to  maintain  the  same  resupply  time  regardless 
of  demands  on  the  system.  An  advantage  of  this  approach  is  the 
ability  to  decouple  maintenance  policies  such  as  service  scheduling, 
number  of  specialists,  and  amount  of  test  equipment,  from  the  stockage 
policies  involving  the  amount  of  spare  stock  provided.  At  the  same 
time,  it  is  difficult  to  identify  the  end  benefits  of  certain  mainte¬ 
nance  policies  when  such  loose  coupling  to  inventory  calculations 
exists.  (Clearly,  there  must  be  a  tradeoff  in  amount  of  spare  stock 
provided  and  the  service  capability  provided.) 

The  second  assumption  is  generally  violated  because  the  scheduling 
and  priority  for  service  are  typically  based  on  the  current  set  of 
inventory  shortages  (or  at  least  those  shortages  affecting  system 
performance).  Furthermore,  when  system  demands  for  service  are 
affected  by  the  past  history  of  demands  (such  as  when  so  many  aircraft 
are  nonoperational  that  flying  demands  are  reduced  and  there  is  a 
consequent  reduction  in  demands  for  service)  then  the  third  assumption 
regarding  Poisson  demands  is  also  violated.  The  latter  assumption 
may  also  be  violated  when  demands  occur  with  a  higher  variance  to 
mean  ratio  than  expressed  by  the  Poisson  process.  The  violation 
of  the  latter  two  assumptions  generally  means  that  the  models  lead  to 
conservative  predictions  of  behavior  in  that  they  assume  more  demands 
than  might  occur  and  they  assume  longer  service  times  than  might  be 
achieved  given  efficient  scheduling  and  expediting  of  certain  types 
of  service.  In  stockage  calculations,  this  will  lead  to  a  higher 
demand  for  spare  stock  and  in  performance  measurement  with  a  given 
level  of  stock  the  models  will  consequently  estimate  performance 
lower  than  might  actually  be  achieved. 

We  have  examined  an  alternative  modeling  approach  which  involves 
direct  integration  of  the  differential  equations  describing  the 
dynamic  queue  while  assuming  only  the  Markov  property  that  the  next 
state  is  determine  by  only  the  transition  probability  from  the  current 
state.  The  advantage  of  this  approach  is  that  the  dynamic  behavior 
of  queues  as  they  respond  to  server  limitation,  service  rate  changes 
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as  a  function  of  both  state  (number  In  the  queue)  and  time,  and  demand 
changes  as  a  function  of  state  and  demand  can  be  considered.  This 
method  Is  also  appealing  relative  to  Monte  Carlo  type  simulations 
In  that  only  one  Integration  Is  required  to  obtain  the  complete 
time  dependent  probability  distribution  of  number  of  Items  In  a  queue 
with  this  approach  while  the  Monte  Carlo  simulation  may  require  a 
large  number  of  trials  to  obtain  the  distribution  with  a  sufficient 
degree  of  confidence.  Although  the  results  of  this  approach  appear 
promising  there  are  certain  drawbacks  still  to  be  overcome.  First, 
all  "pipeline"  or  different  service  processes  must  be  Included  In 
the  set  of  state  equations  because  of  the  loss  of  the  Poisson  properties 
and,  secondly,  the  states  of  all  components  sharing  the  same  set  of 
servers  must  be  considered  at  one  time.  This  can  lead  to  a  potentially 
large  number  of  state  equations  to  be  solved  because  of  geometric 
growth.  Ve  are  currently  examining  methods  of  reducing  the  number 
of  equations  which  must  be  considered  at  one  time. 
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